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'on : 

O ■ Abstract 
o ■ 

I Sterile neutrinos are studied as subdominant contribution to solar neutrino physics. The 

Qh' mirror-matter neutrinos are considered as sterile neutrinos. We use the symmetric mirror 

■ model with gravitational communication between mirror and visible sectors. This com- 

. munication term provides mixing between visible and mirror neutrinos with the basic scale 

^ I /J, = v^y^/Mpi = 5 X 10~^ eV, where v^w = 174 GeV is the vacuum expectation value of 

^ I the standard electroweak group and Mpi is the Planckian mass. It is demonstrated that 

^ ■ each mass eigenstate of active neutrinos splits into two states separated by small Am^. 

^ . Unsuppressed oscillations between active and sterile neutrinos {ua ^ i^s) occur only in 

O ! transitions between each of these close pairs ("windows"). These oscillations are charac- 

^ I terized by very small Am? and can suppress the flux and distort spectrum of pp-neutrinos 

O ' in detectable way. The other observable effect is anomalous seasonal variation of neutrino 

^ ■ flux, which appears in LMA solution. The considered subdominant neutrino oscillations 
z^a ^ can reveal itself as big effects in observations of supernova neutrinos and high 

Oh! energy (HE) neutrinos. In the case of HE neutrinos they can provide a very large diffuse 

^ I flux of active neutrinos unconstrained by the e-m cascade upper limit. 

^ '. 1 Introduction 

P.: 

1.1 Why light sterile neutrinos? 

A sterile neutrino, z/^, is a neutral, spin 1/2 particle which is a singlet under SU(3)x 
SU(2)xU(l) (SM model) group. In supersymmetry, grand unified theories (GUT) and 
superstring models such singlets appear naturally. It is more difficult however to have them 
light. Sterile neutrino always interacts, sometimes very weakly, with ordinary neutrinos. 
This interaction naturally results in small neutrino mixing. It is more difficult to build a 
model with large mixing between sterile and ordinary neutrinos. 

However, there are many models of light sterile neutrinos, which effectively mix with the 
ordinary ones. It is natural to consider the models in which ordinary and sterile neutrinos 
mix, even if this effect is weak. In our view, an investigation of sterile neutrinos should not 
necessarily start from this or that observational motivation; it should rather select a well 
definite theoretical framework and study its implications, with the hope that some of them 
are observable. This is the strategy of our work. 
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It must be said that at present there are not many indications for effective oscillations 
between sterile and active neutrinos. Two flavor oscillation z/g — z/^ is excluded as solu- 
tion to the solar neutrino problem, especially after SNO data p. Oscillation z/^ — > is 
disfavored as solution of atmospheric neutrino anomaly j^I- In recent works |S1 IH El, 
stringent upper limits on the contribution of sterile neutrinos to the solar or atmospheric 
neutrino experiments are obtained. An indirect hint for a sterile neutrino comes only from 
interpretation of the above-mentioned data when combined with those of LSND [T (for the 
difficulties in such interpretation see jH], and P for further discussion). 

However, a small mixing of sterile neutrinos with active ones and/or very small Am? 
remain a viable possibility. This possibility can reveal itself in subdominant processes, or in 
some new phenomena such as neutrino oscillations at very large distances: neutrinos from 
supernovae and high energy (HE) neutrinos from the mirror matter. 

The models for sterile neutrinos include those based on super symmetric theory, with 
or without R-parity violation e.g. piU ITT], on the GUT and string models, e.g. component 
of 27-plet in Eg model |12j or SM singlet with additional U(l) charge [TSj, on the hidden 
sector, modulinos in supergravity ^3], and many other models. 

We believe that any theoretical model for sterile neutrinos should explain why the mass 
of a sterile neutrino is close to that of the ordinary neutrinos. Actually, this is the reason 
why we select the mirror models JHl El QHI for sterile neutrinos: neutrino masses arise 
only from operators with dimensions 5 or larger, and therefore are suppressed by inverse 
power of large masses, exactly as for ordinary neutrinos. 

In this paper we discuss the origin of the mixing of ordinary and mirror neutrinos 
('communication' term) and estimate its impact on oscillations. The most important effect 
is the splitting of the unperturbed mass eigenvalues of active neutrinos to two close-set states 
(see Fig.^. The transition between these states results in the oscillation to sterile neutrinos 
with long wavelength. We consider the applications of our model to solar, supernova and 
high energy neutrinos. 

1.2 Models with mirror matter 

Mirror matter was first suggested by Lee and Yang [fl^ in 1956, who proposed that the 
transformation in the particle space which corresponds to space inversion x —x should 
not be the usual parity transformation P, but P x R, where R transforms a particle (Lee 
and Yang considered the proton) into a refiected state in the mirror particle space. This 
concept was further developed by Salam [20], but in fact this idea was clearly formulated 
only later, in 1966, by Kobzarev, Okun and Pomeranchuk [21jj. In this work it has been 
proposed that mirror and ordinary matter may communicate only gravitationally, and that 
the objects from mirror matter (stars and planets) can be present in the universe. Okun [22] 
considered also the communication due to new very weak long-range forces and discussed 
this interaction for celestial bodies from mirror matter. Since that time mirror matter 
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Figure 1: The double degeneracy between neutrino mass eigenstates of ordinary and mirror 
world {ui and j/^', respectively, with z = 1, 2, 3) is lifted when the small mixing (communica- 
tion) terms are included. The new mass eigenstates, denoted as and z/~, are maximal 
superpositions of z/j and u'f uf = {ui + u[)/\/2 and = {v[ — z/j)/v^. 

has found interesting phenomenological applications and development [23 1211 • been 
boosted in 1980s by superstring theories with x E'^ symmetry. The particle content and 
symmetry of interactions in each of the E^ groups are identical, and thus the mirror world 
has naturally emerged. 

One can describe the motivation for the mirror sector in the following way. 

The Hilbert space of the particles is assumed to be a representation of the extended 
Poincare group, i.e. one which includes the space coordinate inversion x —x. Since the 
coordinate operations, inversion and time shift, commute, the corresponding operations 
in the particle space, and Hamiltonian H, must commute, too. It implies that parity 
defined as eigenvalue of operator must be integral of motion for a closed system. For 
this it can be suggested that 1^ = P x R, where P is the usual operator of reflection and R 
is the operator which transforms an ordinary particle to the mirror particle. Thus parity 
is conserved in the total Hilbert space of ordinary and mirror particles. The assumption of 
Landau [23] was R = C, i.e. we can say that he suggested to use antiparticles as the mirror 
space, but then CP must be conserved which we know is not the case today. 

Mirror neutrinos as sterile neutrinos and various applications of mirror matter have been 
intensively studied during the past several years in the context of explanation of atmospheric 
and solar neutrino problems IT3| 1^ . cosmological problems, including inflation and 
nucleosynthesis |2Sl 123 12HI 1211 HE] , dark matter and galaxy formation [2011231201121], extra 
dimensions [3^ and high energy neutrinos [TH] . 

Mirror matter scenarios have two basic versions. The symmetric version was suggested in 
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the early works and more recently advocated in [SHIE] ■ The Lagrangian which describes the 
particles and their interactions in the visible and mirror sectors, £vis and ^mirr, are perfectly 
symmetric and transforms into each other when x — >■ —x, accompanying by all left states 
transforming into right and vice versa: ^j? and where primes denote the 

mirror states. The vacuum expectation values (VEV's) of the Higgs fields are also identical 
in both sectors. Parity is conserved in the enlarged space of ordinary and mirror states. 
The masses and mixings of neutrinos in the symmetric mirror model is studied by Foot 
and Volkas in ^3] • The two sectors (ordinary and mirror) communicate through the Higgs 
potential and mixing of neutrinos. This is a phenomenological description, in contrast 
to dynamical description of Refs.[Tni HZj, where the two sectors interact gravitationally 
and neutrino mixing follows from this interaction. Neutrino masses and mixings in each 
sector are induced by the usual see-saw mechanism, and mixing of neutrinos of different 
sectors are postulated as e.g. itivlv'^, where mirrors neutrinos are denoted by primes. As 
demonstrated in ^3] the most general mixing terms compatible with parity conservation 
results in maximal mixing of ordinary and sterile neutrinos. 

The asymmetric version was suggested by Berezhiani and Mohapatra p/Tj. They assume 
that while all coupling constants in the two sectors are identical, the VEV's are different 
and break the parity. The ratio C, = v'/v of electroweak VEV's {{H) = v and {H') = v') 
gives thus the scaling factor for ratios of masses in the ordinary and mirror worlds, such 
as masses of gauge bosons, leptons and quarks. The basic communication between the 
two worlds is gravitational. It is taken in the form of universal dimension 5 operators, 
suppressed by the Planckian mass Mpi. Operating inside each world and between them, 
these terms give neutrino masses and mixings. However, to describe the desired neutrino 
masses, the authors assume also additional communication through the singlet superheavy 
fields, which results in the similar dimension 5 operators suppressed by superheavy mass 
A < Mpi. 

A similar model-with asymmetric hidden sector-was studied in Ref.pH]- The commu- 
nication of the two sectors is described by a dimension 5 operator with superheavy mass 
A in the denominator. The neutrinos in this model are found to be maximally mixed and 
mass degenerate. The neutrino masses and mixings are obtained with help of dimension 
5 operators with one scale A, with two different electroweak VEV's, v and v', in the visi- 
ble and mirror sectors, respectively, and using VEV's of two SU(2) singlets ($) = V and 

($') = v. 

1.3 The nucleosynthesis constraints 

There are several dangers to be watched for in the models with mirror matter. They are 
connected with cosmological (big-bang) nucleosynthesis. 

In mirror models the number of massless and light particles is doubled, and this case 
is excluded by cosmological nucleosynthesis if the temperatures of mirror and ordinary 
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matter is the same. A natural solution might be given by making the temperature of mirror 
matter, T', lower than that of visible matter T. It is rather difficult to implement this in 
the symmetric models [121 • asymmetric models ^7] a natural solution is given by the 
different couplings of the inflaton to the visible and mirror particles [2^1 • The inflaton decays 
with different rates to visible and mirror matter, producing thus the different temperatures 
of the two sectors. In the symmetric models the different temperatures can be obtained in 
a two-inflaton model PH]- In this model there are two inflatons, ip and (p' , with identical 
couplings to the visible and mirror matter, respectively. The roll of the inflatons towards the 
minimum of the potential is not synchronized, and the particles produced by the inflaton 
which reaches the minimum earlier will be diluted by the inflation driven by the second 
inflaton. By definition, the first infiaton is the mirror one. 

Having the temperature T' < 0.5T, does not solve the problem of cosmological nucle- 
osynthesis completely. The number of additional effective neutrino fiavors is limited by 
cosmological nucleosynthesis as AA^,^ < 0.2 — 0.3. Even if the initial density of mirror neu- 
trinos is strongly suppressed, they might reappear again with the equilibrium density due 
to oscillation of the visible neutrinos to sterile neutrinos; for a review and references see 
The oscillations might bring the sterile neutrinos in equilibrium with the active ones. 
(Indeed, while u^s oscillate into z/^'s, the missing z/^'s are replenished again by thermal pro- 
duction.) In case of small mixing angles, when Ua and z/^ can be approximately considered 
as mass eigenstates with masses and m^, the non-resonant oscillation Ua Vs occur 
when Am? = ml — > . In this case the limit on Am^ allowed by nucleosynthesis is 
given by [5^ : 



where 6 is the vacuum mixing angle for Va — mixing, and ^ = 3.16 for Ua = J^e, and 
^ = 1.74 for z/(j = u^/r- Using AN^, < 0.2 — 0.3 one obtains from Eq. (jT]) the upper limit on 
Am^ for given mixing angle 6. 

For large (or maximal) mixing, the upper bound on z/^ (where a ^ e) oscillations con- 
tinues to roughly satisfy the scaling law previously given, whereas the bound on Am^^^, 
becomes much more stringent. In a number of calculations [211 EDI E] bounds were obtained 
in the range 



where the effect of Ua Ve oscillations are neglected. 

The model for neutrino masses presented in this paper has oscillations into sterile neu- 
trinos with small Aw? ^ that satisfy the bounds cited above (see Section 1?!^ . However, we 
note that, a 'priori^ these bounds are not unavoidable. As was first remarked in Ref . |35j . 
these limits become much weaker in presence of large lepton asymmetry L = {n,j — ny)/n^. 
For example, in case of z/g Vg oscillation the limit (Q) is replaced by 



Am^ sin^2^ < ^ lQ-^{AN^f eV^ 



(1) 



Ami v< < 10"® - 10"^ eV^ 



(2) 



AmVeV^ <4-102|Le 



(3) 
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The suppression of Ua — >■ z/g oscillation is due to the matter effects, which appear because 
the neutrino potential depends on lepton asymmetry. The scale for large and small lepton 
asymmetry is given by the baryon asymmetry B ~ 10~^°. The lepton asymmetry needed for 
the above-mentioned effect must be larger than ~ 10^'' jSE]- An application for atmospheric 
and solar neutrino oscillations with L ~ 10~^ is considered in The lepton asym- 

metry can be generated by some unspecified mechanism, but to be larger than B ~ 10~^°, 
the lepton asymmetry must be generated after electroweak phase transition, otherwise it 
would be reprocessed by the sphaleron mechanism into too large a baryon asymmetry. A 
mechanism of generation of the lepton asymmetry is neutrino oscillation itself (see for 
the status and references). Lepton asymmetry is generated only in case of small neutrino 
mixing. At maximal mixing, for example, probability of Va oscillation is large, and 

scattering of z/^ provides thus the equilibrium of sterile and active neutrinos. Foot jl21 sug- 
gested the following model where suppression of oscillation is provided by self-generating 
lepton asymmetry. There are four neutrinos z/,-, z/^, z/g and z/^, among them and Vs have 
small mixing, while the other neutrinos are allowed to have the large mixing. Large lepton 
asymmetry L is produced by Vr oscillation, and L suppresses the oscillations of the 

other neutrinos. (The oscillation z/^ u', where u' is the mirror neutrino, is suppressed 
more strongly than z/^ ^ z/^, because of self- interaction of u' j3Z|-) Thus, even in the case 
of large Am^, the z/^ — > Ug oscillation and the resulting nucleosynthesis restrictions are 
suppressed in presence of large lepton asymmetry, either existing or self-produced. 



2.1 Origin of the mass terms 

Our model belongs to the more general framework of the symmetric models of Foot and 
Volkas ^3]. Differently from these authors, we discuss the origin of the mass terms, and 
focus on the models with gravitational communication terms between ordinary and mirror 
neutrinos. 

2.1.1 Gravitational communication terms 

As in the first classical works ^HJ-ES] on mirror matter we shall assume that that mirror 
particles communicate with the visible ones only gravitationally. For the description of this 
interaction we shall use dimension 5 operators |^ I44j . For neutrinos these communication 
term, obtained from the SU{2)l x [/(I) x SU{2y^ x U{1)' scalar, reads: 
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The model 



comm 




{VaL(l))WpR(l)'), 



(4) 
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where Mpi = 1.2 x 10^^ GeV is the Planckian mass, and = e,fi,T, and are the 
neutral component of the electroweak Higgses from visible and mirror sectors, respectively.^ 

After spontaneous electroweak symmetry breaking the Lagrangian Q generates the 
terms, which mix visible and sterile neutrinos, 

Mpi 

where v= 11 A GeV is VEV of the Higgses, which is assumed to be the same for SU{2) and 
SU{2)' groups. 

We assume that the coefficients Xap are of the order of unity, and this is a natural 
assumption, once we consider non-perturbative gravitational interaction as origin of the 
term (^. In this case we have in our model basically only one mass parameter 

^ = t;2/Mpi = 5.0 X 10'*^ eV. (6) 

However, it is important to remark that in some models this parameter might be slightly 
or essentially different: 

1 ) The non-perturbative gravitational mechanism responsible for the communication term 
(0} could have explicit suppression factors, flavor independent or perhaps dependent. For 
instance, in the wormhole model there is an exponential suppression factor exp(— 5), where 
S is the action, which can be approximately expressed through the wormhole throat radius 
i?, as iS ~ MpiR"^. This radius is inversely proportional to Mpi with unknown proportion- 
ality coefficient. The suppression can reach many orders of magnitude 

2) Eq. (0]) can contain some numerical factors such as l/(47r)^, Clebsch-Gordan coefficients, 
etc. 

3) Also, one has to recall that certain parameters of the standard model, such as gauge 
couplings and top Yukawa coupling, are indeed of the order of unity, but many other pa- 
rameters are much smaller; e.g., Vub ~ 3 x 10~^, rris/v = 7 x 10~^, and rrie/v = 3 x 10^^. 
This suggests that some flavor selection rule, or other mechanism, provides explicit suppres- 
sion of Yukawa couplings. A similar (or the same) mechanism could well be acting on the 
couplings A, and result in their suppression or smallness (e.g. note that in the celebrated 
seesaw model couplings similar to A in Eq. (jH) are related to the neutrino Yukawa couplings, 
that are likely to be smaller than unity). 

The most general neutrino mass matrix in the flavor representation can be written as 

(7) 

for flavor states, and the Latin letters 




mass = { t , ^, , \ + h.C. 



^Here and everywhere below we use the Greek letters a, f3, 
k... for the mass states. 
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where 3x3 matrix m is given by £mix from Eq. Q, and 3x3 matrixes M and M' are 
either complex conjugate or identical due to assumed mirror symmetry (see Appendix). 

We shall assume that the content of matrix M in the visible sector is determined by 
interactions inside this sector, e.g. by the see-saw mechanism. When m = the mass matrix 
M in the mass eigenstates basis is M = diag(Mi, M2, M3), where Mj are masses generated 
by the see-saw mechanism. The numerical values of these masses are outside the scope of 
our work. 

All elements of matrix m are approximately equal and are mainly controlled by the 
fundamental scale of our model fi. Due to M' = M (or M' = M*) the diagonalization of 
mass matrix ((Tj) results in maximal mixing of and u'^. 

Before proceeding to detailed calculations let us consider an illustrative example of two 
neutrinos u and u'. The 2x2 mass matrix in this case is given by 



where we assume Mj ^ jj,. When the interaction between the two sectors is switched off, 
/X = and the neutrinos are mass degenerate. With /i taken into account, the mixing is 
maximal sin 26 = 1 and the mass eigenvalues split to mi^2 = Mi ± /i, so that Am^ = -iMiH 
(more precisely, Am^ = 4Re (Mjm*), since neutrino oscillations depend on the product of 
neutrino mass matrix and its hermitian conjugate). The transition between the split levels 
results in z/^ z/^ oscillation with small Am^. 

This feature survives in the three neutrino case, when each mass eigenvalue, Mj, splits 
into two close ones (see Fig. [T]). This provides additional Ua oscillation with small 

Am^ between the split levels. We shall prove that unsuppressed oscillations of active to 
sterile neutrinos exist only between these split states of one level ( "window" ) with the small 
Am^, while the short wave oscillations (large AM^) are suppressed. 

2.2 Mass spectrum and oscillations 

In case of m = we have the active neutrino oscillations described by the matrix M. We 
introduce its decomposition into mass eigenstate matrix Mj as: 



The precise numerical values of masses Mj and of the mixing angles have to be taken from 
the "standard" models of active neutrino oscillations and thus they are outside the scope 
of our work.^ However, for two masses (M3 and M2) there are the lower bounds from the 

^We will consider the case of "normal hierarchy" , which-beside being consistent with all we know-arises 
most commonly in grand unified theories, models with flavor selection rules, etc. See, e.g., |46| . 




(8) 



M = U* diag(Mj) 



(9) 
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data on atmospheric and solar neutrinos, 



M3 > (AM2,^)1/2 ~ 5 X 10-2 eV 
M2 > {AMliY/^ ^ 7 X 10-3 eV 



(10) 



while the lightest mass Mi remains unconstrained, and its value can be (very) small (e.g. 
Ml <^ M2). Also, we are allowed to borrow the mixing angles suggested by data, without 
dwelling on questions of how to justify their precise values. In the three neutrino case, 
it is convenient to rotate away the phases in the masses Mj, since they do not affect the 
oscillations. We will see below that these phases play a more interesting role in our 6 
neutrino case. Note, that with Mj given above Am^ calculated for oscillations z/„ u' 
{a = e, /i, r) in all the three windows of Fig. [H respect the upper limits ^ and 

2.2.1 Oscillations into mirror neutrinos 

Let us come back again to Eq. ((7j), where the matrices M and m are written in the flavor 
representation. But before, we introduce a notation for the expression of the communication 
term in the special basis where the matrixes M = M' are diagonal: 



The general 6x6 neutrino mass matrix can be brought in exact diagonal form using two 
unitary matrices, defined by: 



with i = 1, 2, 3 the indices for mass eigenstates. The expression we get for the mixing 
matrix, that relates the ordinary and mirror neutrinos to the 6 mass eigenstates, uf, is the 
following one: 



This is the master equation, that now we discuss and analyze. If m = 0, we get two 
equalities: (1) U± = U (namely, the two unitary matrices are equal) and (2) M±i = Mj. In 
this case, ordinary and mirror neutrinos are pairwise degenerate but do not oscillate into 
each other as can be verified from Eq. ()13p. Instead, when m 7^ 0, there are deviations from 
both these two equalities that, in turn, produce different types of oscillations. Their origin 
can be traced back to different terms of the matrix m of Eq. (jlip : 

(1) The off-diagonal terms of rh lead at first order to f/+ 7^ f/_ and to no-splitting between 



m = U^{m)U 



(11) 



Ui{M ± m)U± = diag(M±i) 



(12) 




(13) 
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M_|_j and M_j. Using Eq. ()13|1 . we realize that there are oscillations into mirror neutrinos 
connected with the splittings Mf — Mj, i.e. with short wavelength: 



-Pshort(z^a mirror) = ^ J2a' 

EiiUl^iU+a'i - Ul^.iU-o.'i) exp 

(the sum over a' accounts for disappearance into any of the mirror neutrinos). In this ex- 
pression we have differences between ?7+ and ?7_ , therefore the amplitude is linearly small in 
the parameters of m, and this means that the oscillation probability is doubly suppressed. 
For this reason, these effects are in practice negligible. 



U4j 



(2) The diagonal terms of m remove the double degeneracy of Mj-values: 
M±i = Mi±mii (15) 

This leads to oscillations into mirror states with long wavelengths, associated with the scale 
Am^ =4Re(Mim*J (16) 

namely, using again Eq. (fTSj) : 

^iong('^a mirror) = ^ |f/„ip sin^ {~Te~^ '^^^^ 

It should be noted that, when L/E becomes sufficiently large, this expression averages to 
1/2: the effect is large. 

Note that, using the approximation f/+ ^ U- ^ U, which has been motivated above, and 
plugging in formula (fTSjl the equation z/^ = Uaii^i, one immediately obtains that the would- 
be mass eigenstates z/j are actually maximal superpositions of the true mass eigenstates: 

''^ = ^(''t-^n, (18) 

and similarly for z/,': 

^: = ^i^t + ^n, (19) 

With these equations in mind, it is easy to summarize the pattern of oscillations: apart 
from the ordinary flavor oscillations ('short wavelength') we must also take into account 
that the state z/j are not mass eigenstates, but a superposition of maximally mixed mass 
eigenstates, Eq. (fTHj) (see again Fig. (T)). This leads to further oscillations with 'long wave- 
length', associated with the splitting in Eq. (fTH|l . 
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Till now, we have presented formulae for vacuum oscillations (Eqs. IT^ and IT7j) . Oc- 
currence of MSW effect ^7] leads to other cases: for instance, if v^. is adiabatically con- 
verted to z/2, after a sufficient distance, this state will further oscillate in vacuum, since 
V2 ~ {^2 — 1^2) I At sufficiently long wavelengths, this will produce a disappearance of 
1/2 of neutrinos, that is the same situation that was discussed after Eq. (fT7|) . 

We can summarize this subsection with the conclusion that oscillations to sterile neu- 
trinos which remain unsuppressed are caused by transitions between the split mass levels 
with small Amf in one window shown in Fig. ^ 



2.2.2 The scales of oscillation into mirror states 



The mirror model we are considering has 3 new parameters in comparison with the usual 
ones: these are the three 'small' Am^ of Eq. ()16|). Large values of these parameters, 
as compared with AM^^^, are excluded by solar neutrino data. In this respect the most 
dangerous parameter is oc M2. Indeed, due to f/gs ~ the splitting Am| is decoupled 
from solar neutrino oscillations, and Am^ oc Mi can always be made very small because 



of the arbitrary value of Mi. But M2 oc yAM^^^ is generically large and may create the 
problems with the second "window" in Fig. ^ However, there exist the allowed textures 
of matrices m (m) when oscillations in the second window are suppressed, and we shall 
describe below two such examples. 

(1) Let us consider first the very specific case with m = fi diag(l,0,0). We shall 
demonstrate, in fact, that this specific case arises from a class of initial textures of matrix 
m with all elements of order one, as implied by Lagrangian ^ with Aq/j ~ 1. We perform 
rotation of m to m using the usual mixing matrix U with f/gs = sin = 0, with the maximal 
atmospheric neutrino mixing angle, i.e. ip = 45°, and with a large solar angle u, namely 



U 



I 



V 



"p. 



\ 

1 

72 J 



(20) 



where s^j = sinu and c^j = cos a; and the common notation for the angles is u; = 612, 
(J) = Oi3, ip = 623- Using Eq. (fTT|) . we obtain the communication matrix m which has all 
elements C(l), as should be provided by \a/3 ~ 1: 



m 



V 71^2.. 



"75^20; 


78^2.. 


is2 


-is2 




2'^oj 


-is2 


is2 


2 


2'^uj 



\ 



J 



(21) 



This property remains true generically, even for other values of the starting matrix, e.g. 
rh = fi diag(l, i, 3) [i here is a/— 1), and actually, this happens even when rh is non-diagonal. 
In other words, we may have very small or negligible oscillations to the second window. 
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without violating the condition that the elements of m are of order unity. 

(2) As the second example, we consider the case when all the elements of the communication 

term m are exactly equal to one. That is, m has the following texture: 



m = fi 



/111 
1 1 1 



(22) 



V 1 1 1 



Using Eq. ^ and Eq. ^ we get 



m = fi 




(23) 



We have not written the off-diagonal terms of rfi because as explained above they play 
no role in oscillations. Note that for this texture all terms of m are of the same order of 
magnitude. In particular all diagonal terms are of same order of magnitude, i.e. there is no 
strong hierarchy in the diagonal terms as in the first scheme. Therefore the splitting in the 
first two windows is 

Ami ~ 4/icos^co'Mi cos^i 



where AM^^i is the solar neutrino mass squared difference for the large mixing angle (LMA) 
solution, and we have introduced two phase factors and ^2- For typical values of mixing 
angles in the LMA regime, sin^ u ~ 0.3 we get: Am^ ~ 1 x 10^^ eV Mi and Am^ ~ 
5 X 10~^ eV^. Since Mi is unconstrained, we can always choose a sufficiently small value 
for it, so that the splitting of the first level, or Amf, becomes small and irrelevant for solar 
neutrino oscillations. If taken at face value, the splitting in the second level is too large, 
and the net result is too much suppression for solar neutrinos. However, as mentioned in 
Section r2. 1. II there are several ways to escape this conclusion, for instance, with the help of 
suppression of scale n by exp(— 5) from wormholes effects. Or even, it is possible that M2 
and m22, when regarded as complex numbers, lead to a small Am"^ just because of their 
relative phase (namely, the phase ^2 can be close to 7r/2). Note that in this second case, the 
smallness of Am^ can be attributed to the phase of M2, instead than to a specific property 
of the matrix m. 

We conclude that mirror models with gravitational communication in many cases have 
a mass matrix m (or rh) which produces too large oscillation effects, but there are also 
models compatible with present experimental data. 





(24) 
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0.75 




Figure 2: How the LMA survival probability is modified by the the oscillation into mirror 
states. The values of Ami indicated at the curves, in units of 10~^^ eV^. Note the 
sizeable spectral distortion at low energies. 

3 Applications 

In this section, we will consider the effect of oscillations into mirror neutrinos on solar, su- 
pernova, and high energy neutrinos. An important remark is that the atmospheric neutrino 
phenomenology is not affected by the splitting of the third level z/3, since for the values of 
the L/E relevant to the atmospheric neutrino problem, there are no oscillations into mirror 
neutrinos. The atmospheric neutrino problem remains pure z/^ ^ Ur oscillations to a very 
good approximation, since CHOOZ |3H| constrains the mixing term f/gs to be small. The 
limit of CHOOZ implies also that solar neutrino oscillations are almost decoupled from 
"atmospheric" frequency AM^j^ in the usual three-flavor context. The same remains true 
in our model: it is sufficient to consider oscillations with frequency AM^^; or smaller. 

3.1 Subdominant solar neutrinos oscillations 

To illustrate the role of oscillations into mirror neutrinos, we will discuss now how the LMA 
solution of the solar neutrino problems is affected by these oscillations. We emphasize here 
that our aim is not a detailed global fit to the solar neutrino data. Instead we want (1) to 
demonstrate that solar neutrino data are consistent with the presence of subdominant 
oscillations into mirror neutrinos in our model and (2) to calculate the effects produced by 
these oscillations, which occur mostly at low energies. 
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3.1.1 CASE 1: Am? ^ 0. 

We consider here the case of sterile oscillations associated with the splitting of the first 
level. The splitting in the second window is assumed to be small as discussed in Sections 
12.2.21 and 12.1.11 The splitting of the third level does not affect solar neutrino oscillations. 

The electron neutrino survival probability with the MSW effect taken into account is 
given by 

00 00 OOOf- /\ 

Pee = COS UJ COS UJm + siu UJ siu UJm — COS UJ COS UJm siu (25) 

Here, 5 is the phase of vacuum oscillations 

(26) 

where L is the distance between production and detection, E is the neutrino energy, and 
Am? is the mass squared for oscillations into mirror neutrinos. The mixing angle at the 
core of the sun 0;^ is: 

sin 2a; 2V2GFPeE 

tan 2a;™ = , with a = 27 

cos 2uj — a M2 — Mf 

Note that the above expression for the survival probability is true as long as the propagation 
in the sun is adiabatic. This happens in the case of the LMA solution, on which we elaborate 
here. 

There are two distinct behaviors of Pee at low and high energies: 
Low energy regime : If the solar neutrino scale is the LMA mass squared difference then 
at low energies Um ~ uj- This implies that Eq. becomes 

Pee ~ ^ ~ 2 ^^^"^ ~ '"^^^ ^ ^^^^ ^' ^'^^^ 

This can be cast in a more transparent manner: 

p^^ = p^^/^^ - cos^ 00 sin^ 6. (29) 

where PeJ^^ is the standard survival probability at low energies of the LMA solution. 
High energy regime: At high energies uJm ~ f • This implies that Eq. ()25jl becomes 

Pee = Sm^UJ. (30) 

This is the standard survival probability at high energies of the LMA solution. Therefore 
the crucial feature of our model is that its predictions coincide with the standard MSW 
solution at high energies but are affected by the subdominant sterile oscillations at low 
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energies: The standard MSW solution is modified at low energies, and most noticeably at 
pp neutrinos energies. This is evident from Fig. Q which shows the survival probabilities 
for some values of Ami superimposed on the usual LMA survival probability. 

Indeed, there is an upper bound on Ami which follows from gallium data. For the 
calculation, we took the fluxes from [HII and the cross sections from and used the 
average gallium rate as obtained by the Gallex/GNO and SAGE experiments (SHj, 70.8±4.4 
SNU. As unperturbed case we choose the best fit LMA solution M| — Mf = 6.2 x 10"^ 
eV^ and tan^cu = 0.4 as obtained in which is consistent with the value found in other 
calculations, see e.g. j^. As illustrated in Fig. |21 the largest value of Ami allowed by the 
3 sigma range is 10~^^ eV^. (It is curious to note that this coincides with the scale of the 
"just-so^" solution I^HESHSHl to the solar neutrino problem.) 

Two remarks are in order: 

(i) As can be seen from Fig. ^ there is a large difference between the profiles of the usual 
LMA solution and our model in the energy range 0.2 to 0.4 MeV. This is precisely the 
energy range of the pp neutrinos. So future experiments like LENS [HZj which will study 
pp neutrinos in real time should see a large distortion in the pp spectrum. By contrast, 
the usual LMA solution does not predict any energy dependent distortion in this energy 
range. When one decreases the subdominant scale, the effect weakens though remaining 
still appreciable for a reasonable range of Ami- 

(ii) The spectral distortion somewhat diminishes when the mixing angle u increases. This 
is simply due to the cos^ u factor in front of the oscillatory term. But even for the largest 
value of the angle allowed by data, 41°, there is an appreciable effect for a range of Ami. 

An important task is to investigate possible signals of seasonal variations in our scheme. 
Even if the LMA solution will be firmly established by future experiments like KamLAND 

, an interesting signature of our model will be the presence of further seasonal variations 
in solar neutrino experiments, as we now discuss. We first rewrite our oscillation probability 
Eq. ()25|) in a form suitable for studying seasonal variations. 

P^^ = p^^/^^ - cos^ uj cos^ Um sin^ S, (31) 

where S is defined in Eq. and 

L{t) = Lo{l-e cos nt) (32) 

In the above equation e = 1.675 x 10~^ and = 1.496 x 10^^ m. Q = ^ where T = 1 
year, and t is the time since the perihelion. Hence we can write the phase of the oscillating 
term as 

5 = Soil- ecosnt) (33) 

where in analogy with previous notation we define 5o = Ami Lo/4:E. This modulation 
implies 

sin 5 = sin [(5o (1 — £^ cos Qt)] = sin 6o — e6o cos 6o cos Qt (34) 
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Figure 3: How the gallium rate changes with the new oscillation scale (the parameter Ami). 
For comparison, the 3 sigma experimental band is shown by dashed lines. 

in writing Eq. (j34p we have used the fact that 5o£cosr2t <ti 1 which is always true for the 
range of Ami ^'^^ considering, and for all neutrino energies of interest. Using Eq. (|34j) 
we get: 

Pee = PeJ"'"^ " COS^ U COS^ slu^ 4 , . 

+e 60 sin 25o cos fit COS^ U COS^ UJm = (-Pee) + -Pee COS ilt 

Here, (-Pee) is the time independent part of Pee, while Pge is the amplitude of the time 
dependent part. 

We first apply this formula to Borexino. For an experiment which measures the mo- 
noenergetic Beryllium line the rate, as a function of time, has the following form 

P' 

R{t) ocl + obo cos (Qt) with obo = , „ , ,^ (36) 

(Pee) + 1/5 

(a factor ^1/5 comes from the neutral current contribution). The amplitude a-Be of the 
time varying factor is a function of the small scale Ami, ^^"^ evaluated at the energy 
of the Beryllium line. We find that even for the maximum value of Ami which is 10~^^ eV^ 
the amplitude is only 0.001. So there is no significant seasonal variations at Borexino. 

However, since the subdominant scale in our model is much smaller than the usual "just- 
so" scale, there may be some seasonal variations at lower energies. With this motivation in 
mind, we analyze seasonal variations for the pp neutrinos. We will consider two experimental 
situations, the first one a la LENS, when the pp flux can be measured separately, the second 
one a la GNO, when there are contributions to the signal also from other neutrino fluxes. 
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Figure 4: Amplitude a of seasonal variation of gallium counting rate ('GNO' case) due to 
oscillations into mirror neutrinos as percentage of the excentricity e, in a function of the 
(average) counting rate. We also show the amplitude app of seasonal variation for events 
induced by pp neutrinos only ('LENS' case). The curves are obtained varying Amf in the 
range below 10~^^ eV^ for fixed LMA parameters. 

In the LENS case, the modulation of the signal is: 

R(t) oc 1 + ttpp cos (Qt) (37) 

where Upp is given by 

_ J%piE)aiE)P:jE 

J%p{EME){P^^)dE- ^''^ 

Here $pp is the pp flux as given in BP2000 pTl, and a is the cross section for neutrino 
absorption on gallium. The integral goes from the threshold of 0.23 MeV to the endpoint 
of the pp spectrum. For the maximum value of Am^ = 10~^^eV^ we obtain that the 
amplitude is slightly more than 1 %, which is smaller than the term e ~ 1.7 % resulting 
from geometrical modulation of the distance in Eq. (j32|) . but is however non-negligible. 
Two conclusive remarks are in order: 

(i) It should be noted that the new modulation is in phase with the geometrical modulation, 
in other terms the rate is modulated by (2e + app) cos fit. 

{ii) As mentioned above, all these calculations are done with a fixed LMA solution, namely 
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Figure 5: How the LMA survival probability (solid curve) is modified by the oscillation into 
mirror states (dotted curve, Am^ = 1 x lO""*^^ eV^). 

the best fit in absence of mirror neutrino oscillations. If we decrease uj to the minimum 
value allowed by the LMA solution, the amplitude increases a bit to 1.23 %. 

Let us consider now the case of an experiment like GNO. When we take into account 
the contribution due to the other neutrino fluxes (which do not receive seasonal variations 
in our model) the expression for a^p Eq. (jHHj) gets modified to 



where Roth is the contribution to the gallium experiment due to the other neutrino sources. 
So the effect will be scaled down in the actual case of GNO (see Fig. 0]). The message is 
that an experiment dedicated to measuring the pp flux only will be better suited to look 
for this kind of effect. 

In concluding this section we again stress the novelty of the phenomena: Even after the 
MSW solution is established, the model we consider can lead to differences in the counting 
rate, spectral distortions, and/or seasonal variations for pp neutrinos. 

3.1.2 CASE 2: Am^ ^ 0. 

Let us assume that Amf is strongly suppressed, but Amg is not, see Sections 12.2.21 and 



! %,{E)a{E)PidE 



(39) 



a = 



J%p{E)a{E) {Pee) dE + Roth 



18 



The electron neutrino survival probability becomes 

00 00 OOOf- /\ 

Pee = COS uj COS uJm + siu uj siu uj^ — siu uj siu ujm siu (40) 

(we use the notation of Eq. (j26p for the phase of oscillation with Am^ replaced with ^rn^). 
In this scheme at high energies when uj„i ^ ^ we get 

Pee = Sin^ UJ cos^ 5 (41) 

So the usual LMA behavior is modified also at high energies. When we select the range 
< 10~^^ eV^ there is no spectral distortion at boron energies, but just a scaling down 
of the usual LMA survival probability. This upper bound is obtained by allowing for a 10 
% decrease in the survival probability. Note that the suppression of boron neutrino flux 
will be accompanied by a disappearance of total neutrino flux, which is given by the factor 
$ = cos^^^'', since only half of the z/2 neutrinos reaches the detector. For this reason, the 
ratio of charged- current to neutral-current events yields the same uo that will be obtained 
by terrestrial experiments, say, by KamLAND. For the same reason, one can use the boron 
flux to extract information on 5 only if the absolute flux value is known from the theory, 
and this makes the investigation difficult. 

In this scheme, the gallium rate is weakly dependent on the new scale unlike the previous 
case. In fact, when Amg changes from 10~^^eV^ to 10~^^eV^ the gallium rate diminishes 
only by 3 to 7 SNU. We plot in Fig. the modification at low energies of the usual LMA 
survival probability due to oscillation to mirror states. One sees that even for the largest 
value of Am\ allowed the effect is weak, unlike the case of the first scheme. This is simply 
due to the fact that for typical mixing angles in the LMA region sm^u <C cos^tu. 

The analysis of seasonal variations is exactly analogous to the analysis in the first 
scheme, except that we now use Eq. fl40|) for the analysis. We find that the coefficient a of 
modulation varies between ±25 % in the allowed Am| range. In conclusion, this scheme 
offers much less clear "smoking- gun" signals and its experimental investigation is difficult. 



3.2 Supernova neutrinos 



The framework outlined above has a natural application to supernova neutrinos. In fact, 
the condition that the phase of vacuum oscillation is large, Am?L/4:E ^ 1 can be rewritten 

as 



Am^ > 1.3 X lO^^V 



1 kpc 



E 



20 MeV 



(42) 



where Am? is the mass scale that gives rise to oscillations into mirror neutrinos. We recall 
that the energy of supernova neutrino events lies certainly in the range 1 < E < 100 MeV 
(lowest bound being mostly due to detector characteristics), and the distance of the galactic 
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center is Lg.c. ~ 8 kpc. We will consider two cases where mixing with mirror neutrinos leads 
to observable consequences. The first case is based on rather standard astrophysics of core 
collapse supernovae and will be discussed at length; the second case starts from a bolder 
speculation, about the existence of mirror supernovae, however it has a nicer signature. In 
a sense, these two cases correspond to the classification of oscillations into "disappearance" 
and "appearance" that we follow below. 

3.2.1 Disappearance of supernova neutrinos 

Let us consider neutrinos from a core collapse supernova. Various patterns of oscillations 
into mirror neutrinos are possible, according to which vacuum oscillation develops (i.e. 
which Am? is sufficiently large). Here, we will focus on the simplest possibility: all three 
Amf satisfy the inequality given above and therefore, averaged oscillations into mirror 
states take place for all states. As shown above, see Eq. (fT7|) and following discussion, 
this situation leads to the very simple result that half of active neutrinos of any type Ve,^l,T 
and their antineutrinos reach the detector. Hence the signature of this scenario is that the 
energy observed is half the energy emitted. 

The theoretical uncertainties are the key issues to verify or contradict the predictions of 
the model we propose. It is unclear whether the theoretical predictions of the energy emitted 
in the gravitational collapse Eth will become accurate enough in future to reveal a difference 
in neutrino energy by a factor of two, say: Eth = 4 x 10^^ erg versus Eobs = 2 x 10^^ erg 
(where we assume that oscillations into mirror neutrinos do take place). Indeed, within the 
existing theoretical uncertainties [^01; these two values are only marginally distinguishable, 
and can be found by varying the parameters of the equation of state of the nuclear matter 
within the existing uncertainties. 

For what regards Eobs, we recall that a future galactic supernova exploding at a distance 
~ Lg.c. will yield several thousands of neutrino events at the Super-Kamiokande detector; 
similarly, very significant statistics will be collected at LVD, Baksan, SNO, etc. This should 
permit a determination of the total energy Eobs at much better than 10 %; possibly, without 
assuming energy equipartition of the various neutrino fiuxes, but testing it by future data 
(see below). 

In this respect, the ~ 20 neutrino events collected from SN1987A might seem already a 
statistically useful indication; let us look at the point more closely. In a recent work [UT] , 
Loredo and Lamb find as optimal fit values Eobs ~ 3 x 10^^ erg in practically all models 
with an accretion and a cooling components of the neutrino signal, which are expected in 
the 'delayed scenario' for supernova explosion. On examination of their Fig. 10, one gets 
convinced that oscillations into mirror neutrinos would not contradict present theoretical 
expectations, even though this would suggest a stiff equation of state for nuclear matter. 
Note incidentally that this paper, though being one of the the most thorough existing 
analyses, assumes strict equipartition, does not include oscillations into active neutrinos. 
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and considers only Vf,p — >■ ne"*" signal using the 'leading order' cross section. Equipartition 
is particularly important: for example if Ve carry 1/4 rather than 1/6 of the total energy, 
which is a reasonable value [U^) the best estimate of ^] reduces to Eobs ~ 2 x 10^^ erg. 

3.2.2 Appearance of supernova neutrinos 

Mirror matter and mirror stars are expected to exist \2\\ 123 HZ] , and this could ultimately 
result into an explosion of a mirror galactic supernova. The mirror neutrinos will oscillate 
into the active one, and half of the the original flux will become observable. All other 
radiation of mirror supernova is undetectable, that gives the basic signature of the event: 
there will be no optical burst in the direction of neutrino burst, and no radio or infrared 
radiation will be detected simultaneously with neutrino burst and later. In the case of an 
asymmetric gravitational collapse, another detectable signal is the gravitational radiation, 
see e.g. jUni- An additional signature can be given by high energy mirror neutrino radiation 
from the young supernova shell [B^. 

3.2.3 Remarks on neutrino spectra and "equipartition" 

On the top of the new effects outlined above, and in both cases considered in Sects l^^TTl 
and I3.2.!4 we will have also the (usual) effects related to flavor oscillations. To be specific, 
the fluxes of Vei and active neutrinos at the detector are: 



where F^^ Fg, Fj^ = = = = F^ are the fluxes without oscillations. For defi- 
niteness, we assumed a normal mass hierarchy of neutrinos, and an angle > 1°: These 
conditions leads to adiabatic MSW conversion of {ue) into jz/a) and of {I'e) into (terms 
order 0^ < few % are neglected). These effects-but without the oscillation into mirror 
neutrinos-have been discussed previously by a number of authors jHSl EHl EZI EHl EHI , with 
the generic conclusion that flavor mixing might lead to 'hotter' z/g and z/g fluxes due to 
the F° component.'^ The flux of active neutrinos F/vc^from neutral current events-is un- 
changed; a peculiar signature of mirror oscillations is that even this flux is reduced by a 
factor of 2. 

Note that, by the time when the next galactic supernova will explode, we might already 
know the usual parameters of oscillations, and therefore the usual type of supernova neu- 
trino oscillations. In this connection, we believe that it is important to stress an important 

•^This statement is intentionally vague, in view of the different conclusions reached by those who inves- 
tigated SN1987A neutrinos including oscillations, see e.g. [7D], and in view of the fact that a full fledged 
theory of supernova explosions is still missing I72L I73L 1^ . 
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consideration (important, even if mirror neutrino oscillations do not take place): Using the 
large i/g-produced data samples along with the presumably smaller samples produced by 
neutral currents and Ve-, "we will have chance to test the hypothesis of equipartition by the 
experimental data, or in other words, we will be able to reconstruct the total energy Eobs 
with a minimum theoretical bias. 

Mirror neutrino oscillations for different energy and situations (supernova remnant) has 
been recently discussed in Ref. [74J, stressing in particular the possibility of a distortion of 
the spectrum. Of course a similar possibility exists for the neutrino spectrum from core 
collapse, for the particular case when the parameters of oscillations into mirror states satisfy 
the relation AmfL/E ~ 1. 

3.3 High energy neutrinos 

High energy (HE) neutrino astronomy includes the wide range of energies from E ~ 
100 GeV up to ~ 10^'^ GeV. We will consider here the diffuse neutrino fluxes at very 
high energies. These fluxes can be produced by three principal sources: accelerator sources, 
topological defects (TD) and superheavy relic particles. 

There is a very general cascade upper limit on HE diffuse neutrino flux [7^1 ITU] . It is 
based on the e-m cascade which develops due to collisions of cascade electrons and photons 
with the target photons, e.g. CMB. The cascade is initiated by high energy electrons or 
photons which always accompany the production of a HE neutrino. The diffuse neutrino 
flux /j/(-E) is limited as 

E^UE) < ^uj^as, (44) 

where ujcas is the energy density of the cascade photons left in intergalactic space. These 
photons have the energies in the range of EGRET observations, which give the upper limit 
on the cascade energy density Ucas < 2 x 10~^ eV/cm^. The limit is very general: it 
is valid for all processes of neutrino generation in extragalactic space (e.g. generation by 
TD and by decay of superheavy relic particles) and in the galaxies if they are transparent 
to gamma radiation. 

The only class of sources that escape the cascade upper bound (HH) is comprised by the 
so called "hidden sources" An example of a powerful hidden source of HE neutrinos is 
given by mirror matter. As demonstrated in Ref. |^ in some models the density of topo- 
logical defects in the mirror matter can be much higher than in ordinary one. Superheavy 
particles produced by mirror-sector TD and the products of their decays are sterile in the 
visible world, but mirror neutrinos can oscillate into the visible ones. The flux of these 
neutrinos can be higher than what the limit (j44j) allows. 
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3.3.1 Z-bursts 



Z-burst is a beautiful idea [ZZj of generation of Ultra High Energy Cosmic Rays (UHECR) 
through the resonant production of Z-bosons in the collisions of UHE neutrinos with Dark 
Matter (DM) neutrinos, v + z/qm — ^ 2'° — > all. The resonant energy of UHE neutrino 
is Eq = m\/2mi, = 4.2 x I0^'^m~y GeV, where is the mass of DM neutrino in eV. 
Following Ref . ^H] we calculate the number of Z-bosons produced per unit volume and unit 
time, integrating over energy the diffuse neutrino flux F^.{E) coupled to the Breit-Wigner 
cross-section (j{E): 

fiz = 4™,^ J hXE)a{E)dE = 47in,atEoh{Eo) , (45) 

i 

where lu = J2 , = 'oQ cm~^ is the space density of one flavor DM neutrinos, and 
at = ASirfi^Gp = 1-29 x 10^'^^ cm^ is the effective cross-section with Gp the Fermi constant 
and fi, = 0.019 the relative width of Z-decay to neutrino channel. In Eq. (j45p it is assumed 
that neutrino masses are degenerate, mjy. = m^. Using the limit on the sum of neutrino 
masses < 1-8 eV from the spectrum fluctuations derived from 2dF galaxy survey 

we shall assume in the calculations below the neutrino mass rrii, = 0.3 eV as maximally 
allowed. Formula (j45j) is exact. 

In the case of decaying superheavy particles (TD and superheavy relic particles) UHE 
photons dominate in UHECR signals at energy E > 1 x 10^° eV jTHJ- Their flux can be 
calculated as 

I,{E) = j-hzR,{E)Q,{E), (46) 

where R.y{E) is absorption length of UHE photon and Q^{E) is the number of photons with 
energy E produced (via vr" decays) per one Z-decay. Using the observed UHECR flux at 
E ~ 10^° eV, one can calculate the flux of resonant neutrinos I^iEo) from Eqs. and 
PUj) as 4 X 10~^^ cm~^s~^sr~^eV~^, while the cascade limit ()44j) is 5 orders of magnitude 
lower. 

In Ref. jHUj it was suggested that the cascade limit can be evaded, if X-particles decay 
exclusively to neutrinos, i.e. X — > uu. However, in a recent work jHT] it was demonstrated 
that this decay results in the electroweak cascading in which electrons, photons and pions 
are efficiently produced and thus the cascade limit is valid for this case too. 

3.3.2 High energy neutrinos from oscillations 

The oscillations of mirror neutrinos into the visible ones are characterized by oscillation 
length Lose ~ E/Am^, much shorter than the typical cosmological distance L ~ 100 Mpc. 
The only exceptional case is given by oscillation of the resonant neutrinos with Eq ^ 1 x 
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10^^ GeV in the first "window" (see the Fig.^, where Ami can be as small as 1 x 10 eV^. 
Therefore, the average suppression due to oscillation length is given by factor i. 

The conversion of the sterile neutrinos into visible ones occurs through two stages. Let 
us consider a sterile neutrino u'^ born with a flavor a and energy E. On the short length 
scale Lshort ~ E/AM^, where AM^ = Mf — M| is the mass squared difference of the 
unperturbed states, u'^ oscillates into two other sterile flavors, and we have all three sterile 
neutrinos u'f^ with (3 = e, fi, r. On much longer scale Liong ~ E/Am"^, where Am^ is a scale 
of the window splittings, sterile neutrinos oscillate into visible ones. Taking into account 
that suppression factors due to oscillation length is 1/2, we can calculate the probabilities 
Pi^'i, for conversion of mirror neutrino u'^ into visible neutrino z/^, using Eqs. (jl7|) and (j2Up . 
In particular, for conversion of mirror muon neutrino u'^ we obtain the probabilities 

sin^ 2uj 1 sin^ 2lj 

which depend only on the solar mixing angle u. For conversion of mirror tau neutrino 
u'^ one should replace z/^ by u'^ in Eq. (jlTj) . For completeness we also give the relevant 
probabilities for the mirror electron neutrino z/^ conversion. 

sin^ 2c<j 1 sin^ 2^7 

Note, that as follows from Eq. ()17|) the probability of conversion Pu'^u/s summed over all 
visible neutrinos Uf^ is equal to |. It means that for Z-burst production when all neutrino 
flavors participate in the resonant reaction, the total oscillation suppression Pii^'^ — z/) = |. 

As it was already mentioned, in some cosmological models the density of mirror-sector 
TD can be much larger than in the visible sector. The ratio of neutrino fluxes from mirror 
and visible TD's can reach 2 x 10^ ^H]- In our model this ratio is modified by the proba- 
bilities given by Eq. (j47j) . which are different for different modes of oscillations, but after 
summation over final states of visible neutrinos the ratio remains the same as in Ref.l 



4 Conclusions 

Oscillations between active and sterile neutrinos, being excluded as the main channel of 
oscillations both in the solar and the atmospheric neutrinos, can be interesting suhdominant 
processes. In this paper we have considered mirror neutrinos as the sterile ones. Our 
particular model is the "classical" symmetric model of the mirror matter, when apart from 
L ^ R interchange, mirror and ordinary matter have the same interactions, coupling 
constants and VEV's. The communication between these two sectors is gravitational, it 
is described by a dimension-5 operator, suppressed by Planckian mass (see Eq. (jH)). This 
operator can be additionally suppressed by some factors, most strongly by exp(— 5) due to 
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wormhole transitions. The gravitational interaction between ordinary and mirror particles, 
being too weak for heavy particles of the Standard Model, can reveal itself in case of 
neutrinos. When the communication interaction is switched off, the visible and mirror 
neutrinos have identical "standard" mass spectrum Mj {i = 1, 2, 3), being provided 
e.g. by "internal" see-saw mechanism. The communication term (0)) mixes the visible 
and mirror neutrinos , generating a non-diagonal term m (in the form of 3 x 3 matrix) 
in the neutrino mass matrix (jZj). The mixing of active and sterile neutrinos is maximal 
and each unperturbed mass eigenstate with mass Mj {i = 1, 2, 3) splits into two close 
mass eigenstates with masses M±i = Mi ± fhu and Am? = 4Re(Mjm*j). Thus three split 
levels (windows) arise. Our basic observation is that unsuppressed oscillations between 
active and sterile neutrinos occur only in transitions between the split levels of the same 
window. Thus, these oscillations are the long- wavelength ones, being characterized by large 
Lose ~ E/Arn^. These oscillations are unobservable in atmospheric neutrinos and produce 
subdominant effects in solar neutrinos. 

An interesting case (CASE 1, described in Section [3.11) is given by the standard LMA 
solution for active neutrinos with a perturbation term £comm (or matrix m) taken as in 
Eq. ()21|) of Section 12.2.21 In this case, only the splitting in the first window works: Am"^ 
is small and Arrig is irrelevant for solar neutrino oscillations. In this case, the large energy 
effects are practically the same as predicted by the standard LMA solution, but the flux 
and the spectrum of pp-neutrinos are distorted as shown in Figs. |21 and 121 A rather unusual 
prediction is the anomalous seasonal flux variation shown in Fig. |3] The predicted effects, 
especially distortion of the pp-neutrino spectrum, can be detected by the future LENS 
experiment jFFl . 

The mirror-visible neutrino oscillations (t'mirr ^ J^vis) which are predicted to be observed 
as subdominant effects in solar neutrinos can reveal itself as big effects in observation of 
supernova and HE neutrinos. In the former case (supernovae. Section f3.2.1|) . half of the 
expected neutrinos might actually be missed, though it would be crucial to improve on the 
theoretical expectation for total energy relased in neutrinos to interpret this signal precisely; 
another possibility is finding a neutrino signal without any optical counterpart (Section 
I3.2.2|) . In the latter case (HE neutrinos. Section . a large flux of mirror neutrinos can 
be produced by mirror topological defects. Mirror neutrinos can oscillate into the visible 
ones, while all accompanying mirror particles remain invisible. This allows a large diffuse 
neutrino flux unconstrained by e-m cascades or any other restrictions. The probability of 
conversion of mirror neutrino z/^ into ordinary neutrino z/^j, summed over flavors (3 and 
averaged over distance is 

Y.^K'^P = ^' a, /3 = e, /x, r, (49) 

due to the fact that oscillation length Lose ~ E / Am? is much shorter than L ~ 100 — 
1000 Mpc relevant for production of high energy diffuse neutrino flux. 
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A Neutrino mass matrices and transformation prop- 
erties of ordinary to mirror particles 

We shall consider here the restrictions imposed on the texture of neutrino mass matrices 
by the transformation properties of ordinary to mirror particles. 

Let us study the most general case of three neutrino flavors {a, (3) for ordinary and 
mirror neutrinos z/q,(x) and We shall limit to the case of neutrinos with "Majorana 

masses". For this reason, we adopt the formalism of Majorana spinors, which obey the 
relation 



The neutrino mass part of the Lagrangian is given by: 

1 
~2 



(50) 



(51) 



M'^p 1^{x)Plu'^{x) + {M'^pY V^{x)Pru'^{x) 
rriap T\^{x)PlI''p{x) - m*^ T^{x)PRUp{x) 

where the star denotes complex conjugation and Pl,r are the projection operators. It is 
easy to see that this is the most general case, by recalling that for the Majorana spinors 
^Pl,rX = xPl,r^- The parameters of this Lagrangian (apart from the arbitrary phases of 
the fields) are those of the two symmetric matrices M and M', and of the arbitrary matrix 
m, and we can construct from them the neutrino mass matrix: 



M 



M 



m 




(52) 



In the context of quantum field theory, the mirror symmetry implies that the action 
S = J d'^xC is invariant under left^right transformation of the fields (see Lagrangian (j51|) 
as a particular example). 

Two kinds of left aright transformations of the arbitrary fermionic fields \& and 
generically denoted as 



L,R 



R,L 



(53) 
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can provide this invariance. 
(1) The first one is given by 



^(x) ^ C W{x). (54) 

It is easy to convince oneself that this transformation belongs to the class described by 
Eq. : e.g. the ordinary left current is transformed into a mirror right- current. Applying 
this transformation to the neutrino mass Lagrangian (j5ip . we obtain that the action is 
invariant if the following conditions hold: M^fj = M'^^, rria/s = rripa- Then, the neutrino 
mass matrix reads: 

with the conditions that M and m are symmetric. 

(2) Now let us consider the standard case of mirror transformation 

^ 7° (56) 

Again, this transformation converts left- into right-currents, and belongs to the class de- 
scribed by Eq. fl^ . Applying this to the Lagrangian given by Eq. (jSH), and noting that 
the measure of integration df"x is invariant under inversion, we obtain that the action is 
symmetric if the following conditions hold: Map = (M^^)* and map = (w^/3a)*. Then, the 
neutrino mass matrix reads: 

with the conditions that M is symmetric and m is hermitian. 

In most applications considered in this paper one can take M = M* and m = m*, and 
thus the cases (1) and (2) are identical. However, there can be the cases when the two mass 
matrices in Eqs. ()55p and ()57p lead to different physical situations. Consider, for example, 
the one family case, when the element of the matrix m is a small value of order e and the 
element of M is order 1 (that is the case we have in mind in our study). Let us further 
assume that the first parameter is real and the second is a pure phase factor e*^ (we use a 
very similar notation in Eq. ()24|) ). In the first and second case, we have respectively: 



M. = ( f ) . A.. = ( f ) (5* 

Oscillations are described by the combination AiAi''. Neglecting the terms of order of e 
we obtain: 



2 
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in case (1) and 



M2MI 




(60) 



in case (2). The mixing angle is in both cases maximal, but in the first case Am^ can be 
much smaller than e (when the phase ^ is close to ±7r/2), while in the second case it is 
always of order of e. 
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